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Abstract We describe the Resummino package, a C-\ — \- and 
Fortran program dedicated to precision calculations in the fra- 
mework of gaugino and slepton pair production at hadron col- 
liders. This code allows to calculate transverse-momentum and 
invariant-mass distributions as well as total cross sections by com- 
bining the next-to-leading order predictions obtained by means 
of perturbative QCD with the resummation of the large loga- 
rithmic contributions arising in the small transverse-momentum 
region and close to the production threshold. The results com- 
puted in this way benefit from reduced theoretical uncertainties, 
compared to a pure next-to-leading order approach as currently 
employed in the experimental analyses searching for sleptons and 
gauginos at hadron colliders. This is illustrated by using of Re- 
summino in the context of a typical supersymmetric benchmark 
point dedicated to superpartner searches at the Large Hadron 
Collider. 



1 Introduction 

Over the last forty years, theoretical developments and exper- 
imental discoveries in high-energy physics have led to an ex- 
tremely coherent picture, the so-called Standard Model of par- 
ticle physics. In particular, the recent observation of a neutral 
bosonic particle compatible with a Standard-Model-like Higgs 
boson [1,2] represents an impressive success of this theoretical 
framework. However, many questions, such as the stabilization of 
the mass of a fundamental scalar particle with respect to radia- 
tive corrections, remain unanswered. Consequently, the Standard 
Model is widely accepted as an effective theory implied by a more 
fundamental one. Among the plethora of candidates for this new 
physics theory, weak-scale supersymmetry [3,4] is one of the most 
appealing and studied option. It not only addresses the above- 
mentioned hierarchy problem but also provides a solution for the 
unification of the gauge couplings at high energies and explains 
the presence of dark matter in the Universe. 

Experimental searches, especially at the Large Hadron Col- 
lider (LHC) at CERN, for the supersymmetric partners of the 
Standard Model particles are therefore among the main items 
of the experimental program in high-energy physics. Up to now, 
both the ATLAS and CMS collaborations have mainly focused on 
seeking for hints of squarks and gluino, being strongly produced. 
All results are however compatible with the Standard Model ex- 
pectation [5,6]. As a consequence, limits on the masses of the 
squarks and gluino are pushed to higher and higher scales and the 



experimental attention starts to shift towards the pair production 
of the electroweak slepton, neutralino and chargino eigenstates. 

Investigations at a center-of-mass energy of 8 TeV of the 
trilepton golden signature have already led to bounds of several 
hundreds of GeV on the masses of these particles [7,8]. However, 
in contrast to the strong production channels where estimates 
of signal cross sections rely on precise theoretical predictions at 
the next-to-leading order and next-to-leading logarithmic accu- 
racy [9,10,11,12,13,14], searches for the weak superpartners are 
only based on next-to-leading order computations, suffering from 
larger theoretical uncertainties [15,16,17]. 

For an efficient suppression of the Standard Model back- 
ground and a more precise extraction of the underlying super- 
symmetric mass limits (or measurement of the supersymmetric 
parameters in the case of a discovery) , accurate theoretical calcu- 
lations of signal cross sections and key kinematical distributions 
are imperative. Along these lines, while supersymmetric particle 
pairs are produced with a vanishing transverse momentum (pt) 
at the lowest order in the strong coupling ct s , gluon radiation 
from quark-antiquark initial states and their splitting into quark- 
antiquark pairs at O(ot s ) induce transverse momenta extending 
to larger values. Next-to-leading order corrections have therefore 
to be taken into account. However, the perturbative calculation 
diverges at small py, which indicates the need for soft-gluon re- 
summation to all orders and for a consistent matching with the 
fixed order results. On the same footings, when supersymmet- 
ric particle pairs are produced close to the production threshold, 
soft gluon emission again leads to large logarithmic terms that 
must be resummed and matched to fixed order in order to obtain 
reliable predictions. 

These considerations motivate us to introduce in this work 
the computer code Resummino, implemented in C++ and For- 
tran and that can be downloaded from the website 

http : //www . resummino . org 

It combines a leading order calculation of slepton and gaugino 
pair production at hadron colliders [18,19,20] supplemented by 
the associated supersymmetric QCD next-to-leading order cor- 
rections [16,17], with the resummation of the leading and next- 
to-leading logarithms to all orders in the threshold [16,17,21] and 
small transverse-momentum [21,22,23] regimes or simultaneously 
in both [21,24,25]. 

In Section 2 of this paper, we briefly review the three resum- 
mation formalisms included in the Resummino package, giving 
the explicit form of the implemented formulas, as well as the 



techniques employed to calculate the underlying integrals. Sec- 
tion 3 is more technical and is dedicated to the installation of 
the program and its running. A set of illustrative results are then 
shown in Section 4, after choosing a benchmark point typical for 
supersymmetry searches at the LHC. 



2 Resummation formalisms 



inverse Mellin transform is consequently necessary to get back to 
the physical space, which implies to take care of the singularities 
possibly arising at the level of the iV-space cross section. This is 
achieved by choosing an integration contour inspired by the prin- 
cipal value procedure and minimal prescription [26,27]. On this 
contour, the Mellin variable N is parameterized as a function of 
two parameters C and <j> and one variable y, 



N(y) = C + ye ±ilt > with ye[0,oo[ 



(1) 



2.1 Main features 

In perturbative QCD, the doubly-differential cross section de- 
scribing the production, at hadron colliders, of a gaugino or a 
slepton pair with an invariant mass M and a transverse momen- 
tum p F is calculated by means of the QCD factorization theorem. 
The fixed order result, is in this way obtained by convolving the 
partonic cross section da ab , computed at a given order in the 
strong coupling, with the universal densities f a and /;, of the 
partons a and b carrying the momentum fractions x a and x b of 
the colliding hadrons, 



2 dM 2j2 ^=J2j dx a dx b dz\x a f a (Xa,l4J\ 

T ab 

\x b f b (x b ,[i 2 F )^ \zda ab (z, M 2 ,p\, p%, ix R )\s(r-x a x b z) 



(1) 



In this expression, the unphysical factorization and renormaliza- 
tion scales are respectively denoted by /ip and fi R and we have in- 
troduced the quantity r = M 2 /Sh, Sj, standing for the hadronic 
center-of-mass energy. After performing a Mellin transform with 
respect to the variable r, this cross section can be reexpressed 
as a simple product of parton densities and the partonic cross 
section in the conjugate Mellin iV-space, 
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Under this form, it is possible to effectively resum to all orders 
in a B the large logarithmic terms arising in the critical regions, 
i.e., when the transverse momentum tends towards zero and/or 
when the partonic center-of-mass energy is close to the produc- 
tion threshold. The exact form of the resummed cross sections is 
detailed, for the small transverse momentum, threshold and joint 
regimes in Section 2.2, Section 2.3 and Section 2.4, respectively. 
Although these large logarithms must clearly be resummed in 
the critical regions, the full perturbative computation, only par- 
tially accounted for by resummation, is expected to be reliable 
away from these regions. In order to obtain valid predictions in 
all kinematical regions, the fixed order (o-( f - -)) and resummed 
(cr( res ')) results have then to be consistently combined by sub- 
tracting from their sum their overlap (cr( ex P)), 
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The latter is obtained by expanding the resummation formula 
to the desired accuracy, i.e., at 0(ce s ) in our case, and is thus 
dependent on the employed resummation formalism. The ana- 
lytical expression of o-^ exp ' are therefore given in the relevant 
subsection below. 

While the fixed order result <r' f °' can in general be directly 
computed in physical space, or in other words by directly using 
Eq. (1), the resummed component and the calculation of its ex- 
pansion at a given order in a s requires factorization properties 
that only hold in conjugate spaces so that Eq. (2) is employed. An 



The parameter C is chosen such that, on the one hand, the poles 
in the Mellin moments of the parton densities related to the Regge 
singularity lie to the left of the integration contour and that, on 
the other hand, the Landau pole related to the running of the 
strong coupling constant lies to its right. Furthermore, the phase 
<t> can formally be chosen anywhere in the range [7r/2, 7r[. 



2.2 Threshold resummation at the next-to-leading 
logarithmic accuracy 

In the threshold regime, the slepton or gaugino pair invariant 
mass M 2 is close to the partonic center-of-mass energy or equiv- 
alently, the Mellin variable AT tends towards infinity. In this case, 
refactorization allows us to rewrite the partonic cross section a ab 
obtained after integrating Eq. (2) over the transverse momentum 
Pt into a closed exponential form [28,29,30,31,32,33], 



^(res.) 

a b 



(N,M 2 , l , 2 F ,^ 2 R )=H ab {M 2 ,^ 2 F ,^ 2 R ) 
xexp[g a6 (iV,M 2 ,/4'>/4)] ■ 



(5) 



The hard part of the scattering process, independent of the Mellin 
variable N, is embedded within the perturbatively computable 
function ~H. ab . On different footings, the function Q ab , that can 
also be computed by means of perturbation theory, collects soft 
and collinear parton radiation and absorbs the large logarithmic 
contributions arising at fixed order. 

It has however recently been shown that Eq. (5) can be im- 
proved, at the next-to-leading logarithmic accuracy, by includ- 
ing and resumming subleading terms stemming from universal 
collinear radiation of the initial state partons [34,35,36,37]. This 
improvement procedure leads to a resummed cross section ex- 
plicitly depending on the one-loop approximation of the QCD 
evolution operator E^ b , which drives the behavior of the parton- 
into-parton density functions with the energy and encompasses 
collinear radiation [38]. Consequently, the original resummation 
formula is modified to 

a^i^M 2 , „%,„%) = £ E^){N,M 2 /N 2 ,^ 2 F ) 

a',b' 

x E^(N,M 2 /N 2 ,^ 2 F )H a , b ,(M 2 ,f, 2 R ) ( 6 ) 
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where we have introduced the reduced Mellin variable defined by 
./V = Ne lE . The improved soft and hard functions Q ab and 1hL ab 
can still be perturbatively computed and read, at the next-to- 
leading logarithmic accuracy, 
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The arguments of the leading and next-to-leading contributions 
to the Sudakov form factor Q ab depend, in addition to the reduced 
Mellin variable, on the one-loop coefficient of the QCD beta- 
function Pq . This quantity and the associated two-loop coefficient 
/3i given for further references are defined, in our normalization 
conventions, by 



ft = y Ca - l N f T R 



*=|[ 
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3C F N f 



(8) 



for Nf active quark flavors. In those expressions, the group theory 
invariants of SU(3) are standard and read Ca = 3, Cp = 4/3 
and tr = 1/2. 

The first two coefficients eft,' and g\J of the function Q ab 
allow to resum the leading and next-to-leading logarithmic contri- 
butions yielded by soft and collinear radiation. In the MS renor- 
malization scheme, their functional form is explicitly given by [16, 
17] 

^ )(A) = 2ik( Ail)+ ^ 1> )( 2A+ln(1 - 2A) )' 
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(9) 



L[Ai 1 '+4 1 '][2A+ln(l-2A) + iln 2 (l-2A)] 



W$ 



This shows that the first two terms of the perturbative expan- 
sion of the Sudakov form factor only depend on the A a and B a 
resummation functions respectively describing soft and collinear 
radiation and flavor-conserving collinear radiation. The leading 
terms of their expansion as series in powers of ct B have been cal- 
culated in the past and are given by [29,39,40] 
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(11) 



for the A a and B a functions. 

Turning to the hard parts of the resummed cross section, 
the leading and next-to-leading order contributions to the 'H a {, 
function read [16, 17] 
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(12) 



where the dependence on the renormalization scale [i R of the 
infrared-finite part Ao of the renormalized virtual corrections crYu 
is understood. The expression of the H b component is not gen- 
eral and assumes that the virtual contributions to the next-to- 
leading order cross section are normalized, in 4 — 2e dimensions, 



It should be noted that for supersymmetric processes, these vir- 
tual corrections not only include pure QCD contributions, but 
also supersymmetric diagrams with squarks and gluinos running 
into the loops. 

In order to avoid a double-counting of the logarithmic terms 
when combining the resummed cross section with the fixed-order 
result as given in Eq. (3), it is necessary to expand Eq. (6) at 
0(a B ). The expanded partonic cross section is then given, in 
Mellin space, by [16, 17] 
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the quantities 7^ 6 being the full Mellin moments of the one- loop 
approximation of the Altarelli-Parisi splitting functions in four 
dimensions, 
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(15) 



^ } W = c '[>(^-i) 

1 1 N ~\ 

^ )(Ar) = 2C 4jaFrT) + (N + i) { N + 2) -^- k 

and 7^ ' the field anomalous dimensions corresponding, in axial 
gauge, to the virtual, iV-independent, pieces of j aa ' [41], 
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7? = -j- and Tg = ft ■ 



(16) 



2.3 Transverse- momentum resummation at the 
next-to-leading logarithmic accuracy 

In order to be able to refactorize Eq. (2) without performing the 
integration over the transverse momentum px and hence subse- 
quently resum the large logarithmic contributions arising at small 
PT, it is necessary to apply a Fourier transform to the partonic 
cross section a ao , 



roc , 

CW M 2 ,p!,/4,/4) = / Ab- J Q (bp T ) 
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(13) 



where Jo (y) denotes the zeroth-order Bessel function. This opera- 
tion renders the cross section explicitly dependent on the variable 
b, conjugate to the transverse momentum py, dubbed the impact 
parameter. This additional transformation allows us to correctly 
take into account transverse-momentum conservation so that the 



partonic cross section can be rewritten under a form where soft 
and collinear radiation exponentiates [42,43,44], 

a' ,a,b' ,6 

x E$(N,1/P >fi %) C &a ,(N,l/b 2 ) C bb ,(N,l/b 2 ) ( 18 ) 

xH a „ b „(M 2 ,v R ) exp [£,»(// (M 2 6 2 ,M 2 ,/4)] . 

In this equation, holding at the next-to-leading logarithmic ac- 
curacy, the presence of the one-loop approximation of the QCD 
evolution operators allows for evaluating the parton densities at 
the natural scale of the process 1/b, with b = (fe/2)e 7E . Moreover, 
all the other functions can be calculated perturbatively. 

Although there are some freedoms, corresponding to the choice 
of a resummation scheme, in the way to separate the different con- 
tributions into the various C a ij, Q ab and "}i ab factors [45,46], we 
adopt the most physical option where the Sudakov form factor 
and the C ab function are free from any hard contribution. In this 
case, the Sudakov form factor is written as [22,23] 



G ab (M 2 b 2 , M 2 , &) = g™ ^/3b ln[M 2 P) ) ln[M 2 b 2 ] 



■ 9 %\^ ln[M 2 b], 



(19) 



where the first term in this expansion, 

5 ^ )(A) = 2ik ( ^ 1)+ ^ 1))[A+ln(1 ~ A)] ' (20) 

collects the leading logarithmic contributions, and the second 
term, 

»£>(». M'/A) = ji-Jfli" +B<"} l»(l - A) 
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the next- to-leading pieces. We recall that the relevant coefficients 
of the resummation functions A a and B a have already been in- 
troduced in Eq. (10) and Eq. (11). 

In the 'physical' resummation scheme that we have adopted, 
the hard function Hab is f ree from any logarithmic contribution 
and includes, as for threshold resummation, the finite parts of 
the renormalized virtual contributions Ao defined in Eq. (13). It 
reads, at the next-to-leading order accuracy [22,23], 



H 



(22) 



ab (M 2 ^ 2 R )=^(M 2 )[l + ^-Ao] . 

Finally, the C ab functions are evaluated, still at the next-to-leading 
logarithmic accuracy and in this scheme, as [22,23] 

C ab (N,t, R ) = S ab + ^[^-C a 6 ab -~, { a 1 b ) - i (N)\ , (23) 

where 7„ 6 denotes the 0(e) parts of the Altarelli-Parisi split- 
ting kernels in Mellin space, 
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J& U W = 



After resumming the partonic cross section in the impact 
parameter 6-space, the resummed cross section has to be trans- 
formed back to the physical py-space. This procedure requires 
to pay a particular attention to the singularities present in the 
resummed exponent when A = 1 in Eq. (20) and Eq. (21) that 
are related to the presence of the Landau pole in the perturbative 
running of the strong coupling constant. Following the prescrip- 
tion presented in Ref. [47], the inverse Fourier transform is cal- 
culated after deforming the integration contour of the b-integral 
into the complex plane by defining two integration branches 



b-e ±iv t with te[0,oo[ and <p £}0, -[. 



(25) 



The Bessel function Jrj(y) appearing in Eq. (17) is then replaced 
by the sum of two auxiliary functions hi and h^ that distinguish 
positive and negative phases in the complex 6-plane, 
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-iy sin 



For any choice of the ^-parameter, these two functions are always 
finite and their sum is independent of v. This splitting has the 
advantage that each of the two functions is associated with only 
one single branch of the integration contour of Eq. (25). 

In order to match with the fixed order result, making use of 
Eq. (3), the resummed cross section of Eq. (17), together with 
Eq. (18), we expand these two equations at order 0(a B ), 



(cxp.) / 



S^-'vlEK'^.^hi'w 



(27) 



>(N,M 2 , P 2 T ,p 2 F , l t 2 R )=H^(M 2 ^ 2 ) 

+ %-hH\m 2 ,» 2 ) 

M 2 

Up 
r'F c 

+gE[^° ) ( M Vi)4 1) w 

c 

+ £E[ c -w«c°V 2 ,/4) 

c 

-^nHHM 2 ,, 2 R )[A^ + A^i]j 2 

-^n^(M 2 ^ 2 R )[B^+B^]j. 

We recall that the resummation coefficients appearing in this ex- 
pression have already been introduced in Eq. (10) and Eq. (11) 
and that the Altarelli-Parisi splitting kernels have been presented 
in Eq. (15). Moreover, the first two coefficient of the perturbative 
expansion of the hard function % n (, and those of the function 
C ab are deduced from Eq. (22) and Eq. (23). In addition, all the 
dependence on the transverse momentum has been embedded 
within the integral J defined by 



J ■ 
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(28) 



2.4 Joint resummation at the next-to-leading 
logarithmic accuracy 

In this section, we generalize the results of Section 2.3 so that 
both types of large logarithms arising either in the small px re- 
gion or near threshold are resummed simultaneously. Since these 



logarithms have the same dynamical origin, their joint reorga- 
nization is possible. In this way, they eventually exponentiate 
very similarly to the case of the transverse momentum regime of 
Eq. (18) [47,48,49], 



a' ,a,b' ,b 

xE^(N,M 2 / x 2 ,H 2 F )Caa>(N,M 2 /x 2 )C~ bb ,(N,M 2 /x 2 ) 
xH a „ b „(M 2 ^ 2 R )exp[g a „ b „(M 2 ,N,b,H 2 R )] ■ 



(29) 



In order to ensure a proper refactorization of the cross section, 
a Fourier transform has again been performed, as in Eq. (17). 
Furthermore, we have introduced the function \, defined by 



X = x(N,b) 



TV 



1 + b/N 



(30) 



which interpolates between TV in the threshold region, when TV 3> 
5, and b in the small-py region, when b 3> TV. Even though there 
are several ways to define such an interpolation, the choice of 
Eq. (30) first implies that the leading and next-to-leading loga- 
rithms both in b and TV are correctly reproduced, respectively in 
the limits b — > oo and TV — > oo. Next, it avoids the introduction 
of sizable subleading terms into perturbative expansions in a s 
of the resummed formula of Eq. (29) that are not predicted by 
fixed-order computations. 

While the Hat, and C ab functions have exactly the same form 
as their counterparts in the small transverse-momentum regime 
shown in Eq. (22) and Eq. (23), the Sudakov form factor now 
reads 

G ab (M 2 ,N,b 2 ,p 2 R ) = 3^(1^0 In x) In* 

♦«!?(£*-*"».£). <31) 

where the coefficients of its next-to-leading logarithmic accurate 
expansion are, in the MS-scheme, given by [24,25] 

»£ ) w=^(^ 1) +4 1) )[ 2A+ln(1 - 2A) ] • 
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(32) 



We recall that the resummation coefficients have already been 
shown in Eq. (10) and Eq. (11) and that the one-loop approxima- 
tion of the field anomalous dimensions 7„ are the TV-independent 
parts of the Altarelli-Parisi splitting kernels (in axial gauge) given 
in Eq. (16). 

Double counting implied when combining the resummed re- 
sults presented above, after getting back to the physical py-space 
as shown in Section 2.3, is again removed by subtracting the ex- 
pansion of Eq. (17), together with Eq. (29), at the first order in 
a s . This expansion reads has the same functional form as Eq. 



(27), after replacing the integral of the zeroth-order Bessel func- 
tion J by 

r°° b 

J= 6b- Jo(bpr) ln X - (33) 



3 Installing and running Resummino 
3.1 Requirements and technical details 

In order to use Resummino, several external libraries and header 
files are required and must be installed on the system. 

First, information on the benchmark supersymmetric sce- 
nario under consideration is passed to the program by means 
of files compliant with the Supersymmetry Les Houches Accord 
(SLHA) conventions [50,51]. We have adopted the choice to in- 
ternally handle such files by making use of SLHAea [52], a C-\ — \- 
header-only library dedicated to input, output and manipulation 
of SLHA data. While SLHAea is fully included in Resummino 
and hence does not need to be downloaded by the user, this tool 
relies on some headers of the Boost C-\ — \- libraries [53] that are 
in contrast not provided with Resummino. Therefore, their pres- 
ence on the system is a necessary prerequisite and the two Linux 
packages boost and boost-devel must be available. 

Next, both the fixed order and resummed components of the 
hadronic cross section require the evaluation of parton distribu- 
tion functions, either in the physical rr-space or in the conjugate 
Mellin TV-space. The Resummino package does not come with 
any built-in parton density fit and entirely relies on the exter- 
nal LhaPdf library [54] which must therefore be installed by the 
user. On run-time, the Mellin moment of the parton densities pa- 
rameterization under consideration are obtained by a numerical 
fit performed by means of the Levenberg-Marquardt algorithm 
dedicated to multidimensional fits of non-linear functions [55,56], 
as implemented in the Gnu Scientific C-\ — \- Libraries (GSL). As 
a consequence, both GSL header and library files must be in- 
stalled by the user before being able to run Resummino. 

The knowledge of the parton densities both in x-space and 
TV-space allows to compute all three components of the hadronic 
differential cross section d 2 a/dpxdM associated with the imple- 
mented physics processes, as described in Section 2. Additional 
integration upon the invariant mass M of the gaugino or slep- 
ton pair or upon their transverse momentum px then leads to 
the singly differential cross sections dcr/dpT and dcr/dM, respec- 
tively. Furthermore, integration upon both variables allows to ex- 
tract total production rates. Let us note that in the case of thresh- 
old resummation, the integration upon py has been performed 
analytically and there is no way to access doubly-differential cross 
sections. All these integrations, together with the usual two-body 
and three-body phase space integration relevant for the types of 
computations performed in Resummino, are achieved by means 
of an adaptive multi-dimensional integration technique based on 
the importance sampling of the integration domain [57]. To this 
aim, we again make use of the GSL C-\ — \- libraries provided with 
the adaptive multi-dimensional integration Vegas algorithm [58]. 

Concerning the fixed order partonic cross sections, Resum- 
mino is based on the leading order results of Ref. [18] and Ref. [19, 
20] for slepton-pair and gaugino-pair production, respectively. 
Next-to-leading order corrections including both the QCD and 
supersymmetric QCD virtual diagrams are implemented as given 
in Ref. [16] and Ref. [17], the associated finite pieces of the virtual 
loops being computed by means of the QcdLoop package [59]. 
Since the latter is fully embedded within Resummino and thus 
does not need to be installed by the user. 



3.2 Installation 

We recommend the user to always use the latest stable version of 

Resummino that can be downloaded from the webpage 

http : //www . resummino . org 

Once downloaded, the package consists of a compressed tar file 

(resummino-x.x.x.tar . bz2 where x-x-x stands for the version 

number) that must be unpacked, 

tar xf resummino-x.x.x. tar . bz2 

In the case all the prerequisite dependencies of Resummino are 
present on the system (see Section 3.1), it is then necessary to 
generate a Makefile appropriate to the system configuration. 
This is done by issuing in a shell the commands 

cd resummino-x.x.x 
cmake . [options] 

The cmake program checks, in a first stage, that all the de- 
pendencies mandatory for Resummino are correctly installed. 
It subsequently creates a series of Makefile scripts allowing for 
the compilation of the Resummino source files and their linking 
with the dependencies. Two optional arguments can be passed 
to the cmake script. The first of these is related to the LhaPdf 
libraries. In the case they have not been installed in the directo- 
ries referred to by the environment variable LD_LIBRARY_PATH (or 
DYLD_LIBRARY_PATH for MacOS systems), the LhaPdf installa- 
tion directory must be specified by means of 

-DLHAPDF=/path/to/lhapdf 

This instructs cmake that the LhaPdf libraries are stored in the 
directory /path/to/lhapdf /lib and the header files in the direc- 
tory /path/to/lhapdf /include. Equivalently, these two directo- 
ries can be provided separately through the cmake options 

-DLHAPDF_LIB_DIR=/path/to/lhapdf/lib 
-DLHAPDF_INCLUDE_DIR=/path/to/lhapdf /include 

The second optional argument of the cmake script consists of 
information on the directory where the Resummino executable 
has to be created (/path/to/install in the example below). This 
is specified by including the option 

-DCMAKE_INSTALL_PREFIX=/path/to/install 

when issuing the cmake command. 

The Makefile can eventually be executed in order to generate 
a local release of Resummino 

make 

make install 

which can then be further used for physics applications. 



3.3 Running the code 

Once compiled, Resummino can be immediately run from a shell 
by issuing 

resummino filename 

where filename consists of the path to a file containing the set- 
tings of the calculation to be performed. Three extra modes of 
running can also be employed by adding an optional flag when 
executing the code, 

resummino --lo filename 

resummino --nlo filename 

resummino --parameter-log=params .log filename 



The first two choices above allow to respectively compute leading- 
order and next-to-leading order quantities (without matching to 
a resummation calculation). In contrast, the last of the three 
options leads to the generation of a file denoted by params.log 
that includes all the numerical values of the parameters defining 
the supersymmetric benchmark scenario under consideration. 

Now, we turn to the way to encode the computation informa- 
tion in the input file to be parsed when executing the code. First, 
the definition of the collider is passed to the program by fixing 
the nature of the colliding beams and the hadronic center-of-mass 
energy (to be given in GeV). This is achieved by configuring the 
variables collider_type and center_of _mass_energy, i.e., by in- 
cluding in the input file lines of the form of 

collider_type = proton-proton 
center_of _mass_energy = 8000 

Let us note that in the current version of the program, only 
proton-proton and proton-antiproton collisions are supported, so 
that the variable collider_type can only be set to one of the 
values proton-proton and proton-antiproton. 

Next, the two produced superparticles must be referred to 
by setting the variables particlel and particle2 to the rele- 
vant Particle Data Group (PDG) codes [60]. We recall that since 
Resummino is strictly dedicated to the production of the elec- 
troweak superpartners, only slepton, sneutrino, chargino and neu- 
tralino states are allowed as final state particles. For instance, 
the production of a lightest neutralino (whose the PDG code is 
1000022) in association with a negatively-charged next-to-lightest 
chargino (whose the PDG code is -1000037) is encoded as 

particlel = 1000022 
particle2 = -1000037 

The numerical values of the masses of those particles, together 
with these of all the other supersymmetric model parameters, 
are provided by means of a file compliant with the SLHA con- 
ventions [50,51], as already mentioned in Section 3.1. The path 
to this file is specified as the value of the variable slha, 

slha = slha. in 

where in the example above, slha. in denotes a generic SLHA 
file. 

The last pieces of information to be included in the input 
file define the type of computation to be performed and the nu- 
merical precision to be reached. The variable result allows to 
select the observable to compute by setting its value to total 
(total cross section cr using the threshold-resummation formal- 
ism), pt (transverse- momentum distribution da/dpT using the 
p^-resummation formalism), ptj (transverse- momentum distri- 
bution da/dpT using the joint-resummation formalism) or m (inv- 
ariant-mass distribution dcr/dM using the threshold-resummati- 
on formalism). For the last three possibilities, the numerical value 
of the transverse-momentum and the one of the invariant-mass at 
which the differential cross section must be respectively evaluated 
have to be referred to via the variables pt and M. For instance, 
implementing in the input file 



result = pt 
pt = 50 

leads to the evaluation of da/dpj- for pj- 

result = M 
M = 600 

implies the evaluation of dcr/dM for M 
having instead 



50 GeV, while 



600 GeV. Furthermore, 



result 



total 



defines the computation of the total cross section. The parton 
density sets to be employed for both the leading-order and higher- 
order components of the calculated observable are indicated fol- 
lowing the LhaPdf conventions which are based on an ordering 
according to parton density group names and numbers [54]. For 
instance, the command lines 

pdf_lo = MSTW20081o68cl 

pdfset_lo = 
pdf_nlo = MSTW2008nlo68cl 

pdfset_nlo = 

instruct Resummino to use the best fits (the variables pdf set_lo 
and pdf set_nlo are set to zero) of the leading-order and next- 
to-leading-order fits of the MSTW 2008 parton densities [61] (as 
indicated by the variables pdf_lo and pdf_nlo). Factorization 
and renormalization scales are internally set to the sum of the 
mass of the produced particles, up to additional factors that must 
be specified in the input file via the intuitive variables mu_f and 
mu_r. For the sake of the example, central scale choices where 
scales are fixed to the average mass of the produced particles are 
enforced by 



mu_f 
mu_r 



0.5 

0.5 



Finally, the speed of the computation is driven by two param- 
eters, the numerical precision to be reached and the maximum 
of iterations allowed when Vegas is numerically computing the 
various integrals presented in Section 2. These are related to two 
input variables, precision, which takes a real number as a value, 
and max_iter which refers to an integer number. Hence, including 
in the input file the lines 



precision 
max_iters 



0.001 
3 



allows for three iterations of Vegas and demands a relative pre- 
cision of 0.1 %. 



4 Illustrative examples 

To illustrate the usage of Resummino in a practical case, we per- 
form several calculations in the framework of one representative 
constrained scenario of the Minimal Supersymmetric Standard 
Model (cMSSM). As designing an experimentally non-excluded 
supersymmetric scenario is going beyond the scope of this work, 
we refer to an earlier study performed by the LHC Physics Center 
at CERN together with both the ATLAS and CMS supersymme- 
try working groups [62]. This analysis is based on 1 fb -1 of LHC 
data, electroweak precision observables and flavor physics con- 
straints. Its conclusion consists of the proposal of several reference 
points in the cMSSM parameter space to be used for supersym- 
metric searches and phenomenological investigations. We adopt 
their 31 st scenario, where the ratio of the vacuum expectation 
values of the neutral components of the two Higgs doublets tan /3 
is set to 40 and the Higgs supersymmetric mixing parameter fi 
is taken positive. At the supersymmetry-breaking scale, the uni- 
versal scalar mass mo i s fixed to 400 GeV, the universal gaugino 
mass m 1 / 2 to 600 GeV and the universal trilinear coupling Aq to 
-500 GeV. After renormalization group running down to the elec- 
troweak scale, squarks and gluino are found heavy, with masses of 
about 1.5 TeV, with the exception of the lightest stop and sbot- 
tom states. The large left-right mixing inferred by the important 
negative value of Aq indeed lowers their masses to 940 GeV and 
1100 GeV, respectively. In contrast, all the other superpartners 
(neutralinos, charginos, sleptons and sneutrinos) are lighter and 



Table 1 Total cross sections associated with the production 
of any pair of superpartners among the lightest gauginos and 
sleptons (£ equivalently denotes mass-degenerate selectrons and 
smuons) in the context of the LHC collider running at a center- 
of-mass energy of 8 TeV and for the benchmark scenario 31 of 
Ref. [62]. The results are computed at the leading order (LO) 
and next-to-leading order (NLO) of perturbative QCD and then 
matched to threshold resummation at the next-to-leading loga- 
rithmic accuracy (NLL+NLO). The corresponding scale uncer- 
tainties are also indicated and resummed cross sections smaller 
than 0.05 fb are omitted. 



Final 

state 


LO [fb] 


NLO [fb] 


NLO+NLL [fb] 


x?x? 

x%x% 


0.1245j££* 
0.0875+^° 


0. 1605+|;^ 
0.1065+3'^° 


0.1554+°;^ 
0.1043+°^° 


xtx° 2 
X1X2 


4.3674+f ^ 

1.4986+^ 


4.8750+^° 
1.7333^2.^° 


4.8248+°-^ 

1.7111+^ 


xt'Xi 


2.8874+1;^° 


3.3463t||% 


3.3086^^.3% 


'■r'-r 


0.0749™ 

0.0477+^ 
1.1428™ 


0.0868+1;^ 
0.0543+|;^ 
1.3809+2^ 


0.0854+°;^° 
0.0534+°|^ 



lie in the 250-825 GeV range. In the following, we further restrict 
ourselves to the production of the lightest electroweak superpart- 
ners whose masses are approximately given by 



m x1 


= 250 GeV , 


m x% 


= m_ ± =480 GeV , 

Xi 




m ML 


= 565 GeV , 


mg R 


= m AR = 460 GeV , 


(34) 


"If, 


= 295 GeV , 


rrif, 


= 535 GeV . 





Considering the LHC collider, running at a center-of-mass 
energy of 8 TeV, we focus in Table 1 on the largest total cross 
sections associated with the production of any pair of two of the 
particles under consideration. We indicate, in the second column 
of the table, results at the leading-order of perturbative QCD, 
employing the leading order set of the 2008 MSTW parton densi- 
ties [61]. In the third column, we compute next-to- leading order 
predictions, convolving the partonic cross section with the next- 
to-leading order set of the same parton density fit. Finally, in 
the fourth column, these last results are matched to threshold 
resummation. Although this does not imply a sensible change 
in the cross sections, we emphasize the importance of resumma- 
tion by showing the theoretical uncertainties that are obtained 
after multiplying and dividing the central scale, set to the aver- 
age mass of the produced particles (see Section 3.3), by a factor 
of two. Stabilization of the results can indeed be observed once 
soft and collinear radiation is resummed to all orders in a s . At 
the leading-order accuracy, the evolution of the parton densities 
introduces a dependence on the factorization scale through po- 
tentially large logarithmic terms, which leads to an uncertainty 
of about ±10%. Although this specific source of uncertainties is 
reduced at the next-to-leading order, new 0(a s ) diagrams im- 
ply an additional dependence on the renormalization scale. This 
yields a total scale uncertainty of a few percents. Finally, expo- 
nentiation, which allows to account for the dominant higher-order 
contributions within the Sudakov form factor, permits to render 
scale variations under a very good control, at the percent level. 

In Figure 1 , we present transverse-momentum spectra related 
to the production of a selection of light chargino, neutralino and 
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Fig. 1 Transverse-momentum distributions for X2X1 an d ^1 ^1" 
production at the LHC, running at a center-of-mass energy of 8 
TeV, at order 0(a s ) (dashed) and after matching the results with 
transverse-momentum (dotted) and joint (plain) resummation. 



slepton pairs at the LHC. We first compute the predictions at 
0(a s ) (dashed) and we next match the results with resummation 
in the transverse-momentum (dotted) and joint (plain) regimes. 
While the fixed-order results diverge at small transverse momen- 
tum due to the large logarithmic terms that have to be exponenti- 
ated, their resummation leads to a finite (and physical) behavior 
with a pronounced peak in the region where pt J; 10 GeV. In this 
kinematical range, the asymptotic expansion of the two resumma- 
tion formulas are in good agreement with the 0(a s ) results since 
they are all dominated by the logarithms. Therefore, matching 
with resummation as presented in Eq. (3) allows for the regular- 
ization of the next-to-leading order predictions for small values 
of pt- The same matching effects also imply that the resummed 
predictions are sensibly larger than the fixed order ones when the 
transverse momentum of the superpartner pair lies within the in- 
termediate py-range of 20-60 GeV. Finally, it is also shown that 
calculations using transverse-momentum and joint resummation 
agrees reasonably well with each other on the entire pr-range, 
although based on different Sudakov form factors. 



5 Summary 

In this paper, we have introduced the Resummino package, a 
CH — h and Fortran program dedicated to precision calculations 
for gaugino and slepton pair production at hadron colliders. The 
program allows to compute total cross sections, invariant-mass 
and transverse-momentum distributions at leading order and next- 
to-leading order of perturbative QCD. In addition, the results are 
then matched to a resummation of the large logarithmic terms 
appearing at fixed order according to the transverse-momentum, 
threshold or joint resummation formalism. 

We have illustrated the usage of our code by adopting a typi- 
cal supersymmetric benchmark scenario for superpartner searches 
and performing in this context various computations by means of 
Resummino. In the presented selection of results, we have chosen 
to emphasize the major advantages of making use of resummed 
predictions, i.e., a drastic reduction of the associated scale uncer- 



tainties and a regularization of the transverse-momentum spec- 
trum in the small-py region. 
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